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I estimate the degree of local chemical equilibration for n, p and u mesons 
in ultra-relativistic nuclear collisions. The ir and p mesons remain near 
chemical equilibrium in all cases, while the u meson density is typically 
30-50% higher than the equilibrium value for O+Ag collisions at y/s = 20 
GeV. If chemical reactions are turned off, the u meson density is much 
larger, 2.5-3 times its equilibrium value. Thus, uj mesons may provide the 
most sensitive tests for the degree of chemical equilibration in nuclear 
collisions. 
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There exist proposals to measure the QCD transition temperature |]J 
and the lifetime of the hadronic gas [fj] in ultra-relativistic nuclear 
collisions, using dileptons from the p° — u peak. To determine the 
effectiveness of these and similar proposals, one must know how close the 
behavior of the rapidly expanding resonance gas is to that of a static 
system at fixed temperature and chemical potential. More specifically, one 
needs to know how close the ir, p and to distributions are to their values at 
chemical and thermal equilibrium. 

In this Letter, I estimate the degree of local chemical equilibration that 
occurs in a central (low impact parameter) ultra-relativistic nuclear 
collision. This is done by calculating chemical (number-changing) reactions 
for an expanding resonance gas while assuming local thermal equilibrium 
and entropy conservation, using several sets of possible values for allowed 
four-body strong interaction reaction rates. The degree of equilibration 
obviously depends on the assumed reaction rates. As these rates are not 
known, I use plausible values (around 1 fm 2 -c for exothermic reactions). I 
neglect interactions with nucleons because the meson density is much 
higher than the nucleon density in the central region of the collision even at 
the lowest energies that I consider. However, this approximation should be 
tested in a future paper, as the cross sections for meson-nucleon 
interactions are typically somewhat larger than those for meson-meson 
reactions. I use the conventions that H = c = ks = 1. 

I find that the ir and p mesons are always near chemical equilibrium. 
For the highest energy (largest y/s) collisions, the u mesons are also nearly 
equilibrated, but significant departures from equilibrium are predicted for 
collisions such as O+Ag at \fs = 20 GeV. Observation of these departures 
from equilibrium would provide confirmation of the approach to chemical 
equilibrium in ultra-relativistic nuclear collisions, and would give a measure 
of chemical reaction rates in a hot resonance gas. 

There are eight exothermic chemical reactions involving at most four 
particles that are allowed in a system of 7r, p and u mesons. Two are simple 
decays: p — > tttt and uo — > tttttt, with rates V p^ n7T = 0.77 fm~ 1 -c and 
^^7™ = 0.04 fm _1 -c at T = 0. I assume here that these decay rates do 
not change with temperature; this may be a bad assumption, but it 
simplifies the calculation, and the temperature dependence is largely 
unknown. 

The other six exothermic chemical reactions are ujuj — > pp, uju — > titt, 
poj — ► 7rp, pp — > 7tcj, pp — ► 7T7T and ixoj — > nn. The rate constants for these 
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reactions are all unknown, so I estimate them in several different ways. I 
obtain the rate constants for the endothermic reverse reactions using the 
principle of detailed balance (the equilibrium rates must be equal for the 
reaction and the reverse reaction) 

I m 
i=l i=l 

where n® q is the equilibrium density of species x. 

I assume boost-invariant longitudinal expansion of the matter formed 
in the collision, with no transverse expansion, following the simple model of 
Bjorken ||. In the absence of viscosity and discontinuities, the entropy of 
the system is constant. I use this to obtain the temperature of the system, 
T, as a function of proper time, r, treating the system as an ideal gas of 
massive particles. 

At V- r,2^ti(k) + k 2 /3e t (k) 



2n 2 T 



2^9iJ o k dk ^ {k)IT _ - = const. (2) 



Here A is the cross-sectional area of the region of hot matter, is the 
degeneracy of species i, and ei(k) = (k 2 + m 2 ) 1 ^ 2 where rrii is the mass of 
species i. I take g w = 3 (isospin), g p = 9 (isospin and polarization) and 
guj = 3 (polarization), with = 140 MeV and m p = = 775 MeV. 

In this paper, I make the optimistic assumption that all of the entropy 
is generated at the beginning of the collision. I then estimate the final 
entropy using the facts that most of the final state particles are pions, 
which are nearly massless, and that a massless ideal gas has approximately 
3.6 units of entropy per particle. Thus, 

Ar [°° ; 2 t, e i(k) + k 2 /3ei{k) 
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V- f°° u2ji M k ) + k 2 /3ei{k) 

2^9iJ o k dk e6i(fc)/T _ 1 = 3.6 dN/dy, (3) 



where dN/dy is the rapidity density of particles in the final state (both 
charged and neutral). This probably results in an overestimate of the 
temperature, but it is sufficient for the purposes of this letter. 

Once A and dN/dy are specified, the only parameter needed to describe 
the expansion of the matter is the hadronic transition temperature. I 
assume that, above some temperature T c ~ 150 — 200 MeV, the matter is in 
either a deconfined or a chirally symmetric phase. The transition from this 
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phase is modelled as occurring in thermal and chemical equilibrium; this 
provides a reasonable starting point, as equilibration rates are expected to 
be significantly faster in these more symmetric phases than in normal 
hadronic matter. This starting point is also supported by data from pp 
collisions at \/s = 27.5 GeV ||, which is fit well by an equilibrium 
distribution at T = 150 - 155 MeV. 

Once T c is specified, the initial proper time (immediately after the 
transition to a resonance gas), To, is obtained from eq. (0). At To, all 
densities have thermal and chemical equilibrium values: 



The subsequent evolution of the system is described by the following set of 
coupled equations: 



Here final states X are summed over, and there are two obvious additional 
terms due to the reaction -r-r-k — ► uj. Exothermic reactions typically have 
weaker temperature dependences than endothermic reactions, so I estimate 
the exothermic reaction rates geometrically and derive the endothermic 
reactions using the principle of detailed balance, as described earlier. I 
obtain T as a function of r by solving eq. (0). [In practice, it is most 
convenient to work in reverse, calculating r as a function of T.\ 

In Fig. 1, I show results with chemical reactions turned off (R = 0) for 
O+Ag collisions at y/s = 20 GeV and for Pb+Pb collisions at y/s = 200 
GeV, assuming that T c = 200 MeV. For the lower energy O+Ag collisions, 
the 7r and p meson densities are not very far from their chemical 
equilibrium values, but the lj density quickly rises to 2-3 times its 
equilibrium value. In the higher energy Pb+Pb collisions, even the u> meson 
density remains near its equilibrium value. The deviation from equilibrium 
is smaller when the chemical reactions are turned on, so I do not bother to 
show results for the higher energy collisions with reactions turned on. 

I allow the gas to evolve for 80 fm/c for the O+Ag collisions and 500 
fm/c for the Pb+Pb collisions; however, transverse expansion develops with 
a time scale of about 10 fm/c, causing the densities to decrease rapidly so 
that the interactions freeze out quickly after about 10 fm/c. This long 
expansion is reasonable for models without transverse expansion, as the 
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final temperature is 100 MeV, and the system reaches a temperature of 120 
MeV halfway through the evolution. If I allowed transverse expansion, the 
system would not evolve as long; however, the final volume would be the 
same (because of entropy conservation) and the expansion time would be 
proportionately shorter, so the meson densities would be approximately the 
same. 

The simplest assumption is that all exothermic rate constants are 
approximately the same. I thus take 

Rij^ki = R (6) 

for all exothermic reactions. I show results in Fig. 2 for O+Ag collisions at 
y/s = 20 GeV with T c = 200 MeV, taking R = 1.5 fm 2 -c (a m 15 mb). The 
uj mesons are still noticeably further from equilibrium than the n or p 
mesons. This is probably because the u density is low, so they are not 
easily destroyed in two-body processes, and their decay rate is also small. 
This effect does not change much for T c = 150 MeV. 

One obvious improvement on this first assumption is to assume that the 
total reaction rates for all pairs of particles are approximately the same, 
but that the relative probabilities of the different exothermic final states are 
given by the probability that the initial state has the same quantum 
numbers (apart from energy and momentum) as the final state. For 
example, for the reaction ujuj — > X, the possible final states are coco, pp and 
tttt. The relative probabilities of these final states are 1:1:1/9 respectively, 
as uu and pp can always be produced from uu, while for nn the vector 
polarizations must combine to produce a scalar final state. The exothermic 
reaction rates are then 



R^ pp = 9R tot /19, (7) 

Rlolo — >7T7r = -R t0t /19) (8) 

p = ^ t0 V4, (9) 

R pp ^ = R tot /12, (10) 

Rpp^mr = R t0t /12, (11) 

Rn^nTT = [l-(n^/An^)]R tot /A. (12) 



R nw ^ n is temperature dependent because one of the possible final states is 
pp, and this reaction is endothermic. 

In order to put an upper bound on the size of the expected effect, I use 
R tot = 6 fm 2 -c (a ~ 60 mb - approximately equal to the pp cross-section 
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[0]). I show the results in Fig. 3 for O+Ag collisions at y/s = 20 GeV with 
T c = 200 and 150 MeV. The oj mesons are even further from equilibrium 
than those shown in Fig. 2, but this is not surprising as the reaction rates 
(and thus the equilibration rates) are smaller. 

This provides a signal for the existence of a hot, nearly equilibrated 
resonance gas - the size of the final state enhancement of uj mesons with 
respect to p° mesons (as they are almost degenerate and thus have almost 
the same density in thermal and chemical equilibrium). This signal is large 
enough to be possibly measurable in collisions of small nuclei at high (but 
not too high) energies, and is probably significantly smaller than the value 
for a non-interacting system (with R = 0). It may be difficult to observe, 
however, as measurement depends on being able to separate final-state 
mesons from those that decay during the earlier stages of the collision. The 
most likely way to see them is to reconstruct the mesons from the 
final-state pions, but I will not attempt to consider the feasibility of this 
procedure as it depends on details of the detectors used for the 
measurement. 

One possible criticism of this signal is that I have neglected changes in 
the uj width at finite temperature. However, this is certainly not true, as 
the chemical reaction terms I use are equivalent to collisional broadening. 
The magnitude of the collisional broadening that I use here is a bit larger 
than the calculated change in the width [|7| , so the results of this simulation 
suggest that it should be possible to see the change in the uj width As the uj 
excess exists throughout most of the collision, it may be possible to 
measure the integrated uj and p densities rather than the final-state 
densities. This could be accomplished with the dilepton signals, if the two 
peaks could be separated. However, this is probably at least as difficult as 
detecting the resonances from analysis of final-state pions, so I again leave 
the feasibility of this technique as an open problem for the experimenters. 

The reported enhancement of <ft mesons in nuclear collisions [|8] may be 
related to the departure of uj mesons from equilibrium that is seen here, as 
shown by Heinz and Lee JJ. The <fi lifetime (50 fm/c) is even longer than 
that of the uj, so it will be further from equilibrium if there are no chemical 
reactions. However, I have not investigated systems including mesons, so 
I do not have any quantitative results at the present time. 

Finally, while ir and p mesons are likely to be nearly in chemical 
equilibrium even in collisions of relatively light ions at relatively low 
energies, this is not true for uj mesons. The departure of the uj mesons from 
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chemical equilibrium is small for collisions of large nuclei at high energies, 
but care will be necessary when comparing to current results at relatively 
low energies. However, even under the most pessimistic conditions 
considered (O+Ag collisions at y/s = 20 GeV with R = and T c = 200 
MeV), the n and p meson densities were fairly close to their equilibrium 
values (factors of roughly 1.3 and 1.6 higher respectively). 
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Figure captions 



1. R = 0, T c = 200 MeV: (a) O+Ag at = 20 GeV (r = 4.6 fm/c), 
(b) Pb+Pb at v/i = 200 GeV (r = 27.8 fm/c). 

2. = 1.5 fm 2 -c, T c = 200 MeV: O+Ag at = 20 GeV (r = 4.6 
fm/c). 

3. i? tot = 6 fm 2 -c: (a) O+Ag at y/s = 20 GeV for T c = 200 MeV 
(t = 4.6 fm/c), (b) O+Ag at = 20 GeV for T c = 150 MeV 
(t = 15.8 fm/c). 
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Figure 1b 
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Figure 2 
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Figure 3a 
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